Superconducting proximity effect in semiconductor nanowires 
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We theoretically consider the proximity effect in semiconductor-superconductor hybrid nanos- 
tructures, which are being extensively studied in the context of the ongoing search for non-Abelian 
Majorana fermions in solid state systems. Specifically, we consider the dependence on the thick- 
ness of the semiconductor in the direction normal to the interface, a physical effect that has been 
uncritically neglected in all prior work on the subject. Quite surprisingly, we find the completely 
unanticipated result that increasing the semiconductor thickness leads to a drastic suppression of 
the induced superconducting gap due to proximity-induced interband coupling. As a result, in the 
limit of strong semiconductor-superconductor coupling, the proximity-induced gap becomes much 
smaller than the bulk superconductor gap and depends weakly on the interface transparency. 



Introduction. The concrete prediction [TH5] about the 
existence of zero-energy localized Majorana modes in 
proximity-coupled semiconductor-superconductor (SM- 
SC) hybrid structures has triggered a great deal of 
research activity aimed at realizing the Majorana 
modes and understanding the details of the underlying 
physics [6H3- Very recent experimental reports [T0UT4] 
providing observational evidence for the possible realiza- 
tion of the predicted Majorana zero-mode in semicon- 
ductor nanowires with proximity-induced superconduc- 
tivity [3H5] have only accelerated the research activity 
in this field, leading to hundreds of papers in the last 
2-3 years. In particular, the nature of the experimen- 
tally observed zero-bias peak in the differential tunneling 
conductance [T0UT2] , a theoretically-predicted |3l H5UT7] 
signature of the elusive Majorana mode, is currently ac- 
tively debated in the literature [T8H29] . 

The theoretical work presented in this paper takes a 
step back from the hot issue concerning the nature of 
the observed zero-bias conductance peak (ZBCP) in the 
topological superconducting phase and addresses a more 
basic question: What are the characteristic features of 
the superconducting proximity effect in SM-SC nanos- 
tructures and how does proximity-induced superconduc- 
tivity depend on the geometrical details of the SM com- 
ponent, in particular on the SM thickness in the direction 
perpendicular to the interface? We believe that a deeper 
understanding of the nature of proximity-induced super- 
conductivity itself is essential to making further progress 
in this field. It is surprising that, although many effects 
that depend on details of the SM-SC structure (e.g., var- 
ious types of disorder, multiband occupancy, length of 
the SM wire and overlap between the Majorana modes 
localized at opposite ends, soft gap, thermal broadening 
of the Majorana mode, decay of the zero-mode due to 
coupling to the tunneling contacts, possible braiding ar- 
chitectures, etc.) have been discussed in the literature in 
some depth, the invariably present finite thickness effect 



corresponding to the finite dimension of the SM layer 
in the direction transverse to the SM-SC interface has 
never been accounted for. In this work, we develop a 
theory for the proximity effect in SM-SC hybrid nanos- 
tructures that explicitly takes into account the multiband 
nature of the SM spectrum due to the finite thickness of 
the SM layer, L z = N z a, where a is the lattice constant 
and N z the number of atomic monolayers in the direc- 
tion normal to the interface, and incorporates the effects 
of proximity-induced inter-band coupling. We empha- 
size that earlier theories have considered the unphysical 
limit L z — 0, when the SM has no band structure arising 
from the ^-confinement, while, in reality, the experimen- 
tal SM-SC hybrid structures used so far in the laboratory 
have L z « 50 — lOOnm. 

We find the qualitatively important and unexpected 
result that, for finite L z above a critical value Lo that 
depends on the parameters of the materials (typically 
Lo ~ 40nm for InSb- and InAs-based structures), the 
proximity-induced SC gap becomes strongly suppressed, 
depends non-monotonically on the effective SM-SC tun- 
nel coupling 7, and vanishes in the strong coupling limit. 
This is in complete contrast with the thin layer limit 
L z — >• 0, the only case studied so far in the literature, 
where the proximity-induced gap increases monotoni- 
cally with the effective SM-SC coupling up to the maxi- 
mum value given by the bulk gap of the SC system. The 
mechanism responsible for this surprising result is gener- 
ated by the proximity-induced coupling among the SM 
bands produced by the quantum mechanical confinement 
in the direction transverse to the interface. When L z 
is large enough, so that the confinement-induced inter- 
band gaps become comparable to the effective SM-SC 
coupling, the inter-band coupling has a very strong, non- 
perturbative effect on the induced SC gap, and, as a re- 
sult, the induced gap is suppressed. 

The finite transverse thickness-induced suppression of 
the SC gap may provide an explanation for why the 
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proximity-induced gaps corresponding to specific SM/SC 
materials measured in different samples and different lab- 
oratories have almost identical values. This is a rather 
puzzling situation from the perspective of a L z — » the- 
ory: The induced gap is much smaller than the bulk SC 
gap Ao, hence weak effective coupling across the SM-SC 
interface is required, 7 « A . However, in this limit the 
proximity gap is proportional to 7, which is expected to 
show substantial sample to sample variation. Our current 
work shows that in SM systems with Lz > Lq a prox- 
imity gap much smaller that Ao can occur in the strong 
coupling limit, where the dependence of the induced gap 
on 7 is weak, leading to an apparent universality of the 
proximity effect in samples of similar thickness. 

Model and theory. We use a minimal model consist- 
ing of an infinite SM wire along x-direction with rect- 
angular cross section L y x L z proximity coupled across 
a z = interface to an s-wave superconductor that fills 
the z < half-space. Our goal is to explicitly include 
the size quantization effect arising from finite L z into the 
theory of the SC proximity effect. Since our main result 
is independent on L y: we focus on the quasi-2D limit, 
L y — » 00. In the presence of proximity-induced super- 
conductivity, the low-energy physics of the SM nanowire 
is described by a Green function matrix with an inverse 
given by [3 [3D] 



[G x ]w(^) =uj-H n 



(1) 



where H nn ^ with n = (n y ,n z ), is the effective 
low-energy Hamiltonian for the SM nanowire in the 
basis of the transverse wave functions <p nw (w) = 
^2a/L w sm[wn w 7r/L w ], with w G {y,z}. In Eq. (II 



S nn / (u) is the proximity-induced self-energy given by |9 , 
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where o~i and are Pauli matrices in the spin and Nambu 
spaces, respectively, and £ is a constant that depends on 
the details of the SC band structure. The effective SM- 
SC coupling at the interface is represented by the ma- 
trix 7 W = #n z (0)^(0)(5n„n; 5 where g = v F \t\ 2 is a 
constant that depends on the surface density of states 
of the SC in the normal state (z/p) and on the trans- 
parency of the SM-SC interface (t). We have assumed 
that the coupling at the interface is uniform, i.e., inde- 
pendent of position in the ^/-direction, which results in a 
block-diagonal structure of the self-energy with respect 
to (n yi n'y). Note that all confinement-induced modes n z 
perpendicular to the interface are coupled. The eigenval- 
ues of the low-energy states can be obtained by solving 
the Bogoliubov-de Gennes (BdG) equation 



det[G _1 (u;)] =0, 



(3) 



where the frequency is restricted to values inside the 
bulk SC gap, \uo\ < A . To the best of our knowl- 
edge, all existing theoretical and numerical work for 
the SM-SC hybrid nanostructures is based on solu- 
tions of Eq. ([3| that neglect the inter-band coupling 
ln z n' z - While this decoupled band approximation is ex- 
pected to be valid in the strong ^-confinement limit, 
when the confinement-induced inter-band gaps are much 
larger than the proximity-induced inter-band couplings, 
it breaks down for large-enough L z . Note that in the 
decoupled band approximation the matrices describing 
the proximity effect become band-diagonal, e.g., j nn ' = 
Inbnn'- Furthermore, in the static limit corresponding to 
\/a| — uo 2 « |A | in Eq. ([2]), the proximity-induced SC 
pairing potential becomes A^ = jn^o/iln + Ao). This 
limit has been extensively used for discussing the proxim- 
ity effect in SM-SC hybrid systems. We note that, in this 
approximation, A^ « A for j n ^> A (strong coupling) 
and A™ « 7™ for 7™ <C Ao (weak coupling) [30j I3T] . 
In the present work, we relax the strong confinement 
approximation and consider the interband-coupling ef- 
fect arising from the n z bands corresponding to a finite, 
realistic thickness L z . 

The main consequence of relaxing the strong confine- 
ment approximation is that bands with arbitrary n z val- 
ues become coupled, which generates second order pro- 
cesses involving low-energy initial and final n z states 
and high-energy intermediate n z states. The transitions 
between these states are proximity-induced and involve 
hopping into the SC. These processes strongly renormal- 
ize the bare SM band parameters and the induced SC 
pairing potential. In particular, for a given diagonal cou- 
pling strength 7 nn , the inter-band coupling results in the 
collapse of the induced SC gap in thick-enough nanos- 
tructures. Also, we note that, in this regime, the con- 
struction of an effective Hamiltonian that describes the 
SM-SC nanostructure has to involve the high-energy n z 
bands, as they are intrinsically coupled to the low-energy 
bands and renormalize them strongly, which ultimately 
leads to the collapse of the induced SC gap. Thus, as- 
suming that the higher-lying n z bands do not matter, 
as has been universally done in the literature, is simply 
incorrect and leads to qualitatively wrong results. 

Results. To understand the effect of proximity- 
induced coupling involving different transverse confine- 
ment bands, we focus on the matrix structure of the BdG 
equation corresponding to the n z quantum numbers and 
consider a SM thin film - SC heterostructure with uni- 
form coupling across the planar interface. Note that the 
results for the proximity-induced gap will be exactly the 
same for a nanowire with finite width L y , provided the 
SM-SC coupling is uniform. Non-uniform coupling ef- 
fects were previously discussed in the literature [30] in the 
strong ^-confinement approximation. The parameters 
used in the numerical calculations roughly correspond 
to the NbTiN-InSb hybrid system measured in the orig- 
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FIG. 1. (Color online) BdG spectra in the vicinity of k — 
for 2D semiconductor-superconductor heterostructures with 
effective coupling 7 w 0.33Ao and different values of the semi- 
conductor thickness: L z — 20nm (black circles), L z — 40nm 
(orange), and L z = 80nm (blue). Two semiconductor bands 
are partially occupied (n z = 2). The yellow regions corre- 
spond to bulk superconductor states. In the decoupled-band 
approximation, the induced gap A(7, Ao) is independent of 
the semiconductor thickness, but for large L z , proximity- 
induced inter-band coupling determines the collapse of the 
induced gap (blue line). 

inal Delft experiment: effective SM mass m = 0.016m e , 
lattice constant a = 0.65nm, and bulk pairing potential 
Ao = 1.5meV. To limit the number of parameters, we 
consider the case of vanishing spin-orbit coupling and 
zero Zeeman splitting, but we checked numerically that 
the results are generic. The numerical results presented 
in this paper are based on a two-band tight-binding 
model of the SM [9 j, but we have also done similar 
calculations using a more elaborate 8-band Kane-type 
model [9 j. The conclusions presented here are not mod- 
ified qualitatively, although the details of the two mod- 
els differ significantly, e.g., the 8-band model, which has 
many more parameters, is characterized by an effective 
mass that depends on L Z1 and requires different tunnel- 
ing parameters for the conduction and the valence bands. 
We also emphasize that the nontrivial dependence of the 
proximity gap on the semiconductor thickness in the di- 
rection transverse to the interface, the main finding of our 
work, represents a generic feature of the superconduct- 
ing proximity effect that remains valid in the presence of 
spin-orbit coupling and applied magnetic field, e.g., in 
the Majorana-carrying SC topological phase of the SM. 

We present our results corresponding to various val- 
ues of the SM thickness L z and different SM-SC cou- 
pling strengths in Figs. [l}|4j Figure [I] already illustrates 
clearly our key new finding: The proximity-induced SC 
gap depends on the thickness of the SM layer and shows 
a strong suppression for L z = 80nm. We emphasize 
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FIG. 2. (Color online) Dependence of the minimum in- 
duced gap on the strength of the effective semiconductor- 
superconductor coupling. In the weak-coupling limit, 7 — ?► 0, 
the minimum induced gap is independent of L z and is given 
by the induced pair potential Ai n d = jAo/(j + Ao). For 
large couplings, the minimum induced gap depends on the 
semiconductor thickness. In thin semiconductors, A increases 
monotonically and reaches the bulk superconductor value 
Ao in the strong coupling limit, 7 — ► 00. By contrast, in 
thick semiconductors the minimum induced gap depends non- 
monotonically on the effective coupling and vanishes in the 
strong coupling limit. 

that, in the decoupled band approximation, the induced 
SC gap corresponding to a given value of the effective 
coupling r )n z n z = 7 is independent of the system size. 
This approximation holds in the case of strong confine- 
ment, as illustrated by the BdG spectra corresponding 
to L z = 20nm and L z = 40nm, which are practically 
identical, but fails for thicker SM layers. The calculated 
minimum gap as a function of the effective tunnel cou- 
pling across the SM-SC interface for L z = 20, 40, 80nm is 
shown in Fig. [2] The widely used proximity gap formula 
A = 7Ao(7 + Ao) -1 represents a good approximation 
only for very small values of L Z1 or for weak effective 
coupling, but completely fails when the separation of the 
confinement-induced energy levels n z becomes compara- 
ble with the proximity-induced interband coupling j nzn > . 
Note that, for L z > 40nm and strong-enough coupling, 
the induced gap is much smaller that the bulk SC gap 
A and depends weakly on 7. Figure [3] explicitly demon- 
strates our key finding that, for large L z (> 40nm), the 
induced gap is strongly suppressed for intermediate and 
large coupling values and vanishes in the limit 7 — » 00, 
while being described by the "canonical" small L z theory 
only in the weak coupling limit, 7—^0. The canonical 
strong-coupling limit, A = Ao, can only be realized in 
thin SM layers. In addition, as shown in Fig. [2j for 
strong SM-SC coupling and L z > 40nm, the induced 
gap is small and depends weakly on 7. We believe that 
this finding explains why the proximity gaps measured in 
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FIG. 3. (Color online) Normalized minimum induced gap, 
5 = A (7, L z )/A(7, 0), as function of the semiconductor thick- 
ness L z and the effective coupling 7. Thin semiconductors 
(L z < 20nm) are in the decoupled-band regime characterized 
by a normalized induced gap S ~ 1 practically independent on 
7. By contrast, the normalized induced gap in thick semicon- 
ductors varies from S = 1 in the weak-coupling limit (7 — > 0) 
to S « at strong coupling. 

different experimental samples tend to have similar val- 
ues. Fig. [3] can be viewed as a proximity 'phase diagram' 
showing the strong suppression of the proximity gap with 
respect of canonical thin limit value in SM nanostruc- 
tures with L z > 40nm. We note that the crossover value 
Lo ~ 40nm is controlled by the effective mass, which 
determines the energy spacing between the confinement- 
induced bands. 

Finally, in Fig. [4] we show the proximity effect in a 
hybrid system with a SC characterized by a position- 
dependent pairing potential Ao(z). We assume that the 
bulk SC pairing potential dies off toward the interface, 
as expected in a realistic physical situation. Specifically, 
Ao(z) = Ao(oo)(l — e~\ z \/ Xa ), where A is the characteris- 
tic length scale (in units of the lattice constant a) for the 
suppression of the pair potential. The commonly used, 
and rather unrealistic, theoretical assumption is that the 
SC pair potential remains constant right up to the inter- 
face, i.e., A = 0. The position dependence of the bulk 
pairing potential for two values of the healing length A is 
shown in Fig. [4^a). The dependence of the normalized 
induced gap A(7, L Z1 A)/A(7, 0, 0) on the effective SM- 
SC coupling strength for L z = lOnm and L z = 40nm 
is shown in Fig. Qb) and (c), respectively. Surpris- 
ingly, in the weak-coupling limit, 7^0, reducing the 
pair potential near the interface (A > 0) results in an 
enhancement of the induced gap with respect to the uni- 
form case A = 0. For intermediate and strong coupling, 
increasing the healing length A results in a suppression 
of the induced gap that is more effective in thin SM lay- 
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FIG. 4. (Color online) Proximity effect from a bulk supercon- 
ductor with position-dependent pairing potential, (a) De- 
pendence of the bulk superconductor pairing potential on the 
distance from the interface. Ao is suppressed on a length scale 
A and vanishes at the interface, (b) Normalized minimum in- 
duced gap, S = A(7, L z , A)/A(7, 0, 0), for L z = 10 nm and 
different values of the characteristic suppression length A. (c) 
Same as in panel (b) for a system with L z — 40nm. The red 
curves (A = 0) correspond to vertical cuts in the diagram from 
Fig. 4. Note that, in the weak coupling limit, suppressing A 
near the interface results in a slight increase of the induced 
gap. For large couplings, the effect of suppressing Ao near 
the interface is stronger in thin semiconductors. 

ers. Adding this finite healing length effect to the finite 
thickness-induced suppression discussed above results in 
an expansion of the 'universal' regime characterized by 
small proximity-induced gaps nearly independent of the 
effective SM-SC coupling toward lower values of the SM 
thickness. We emphasize that, regardless of L Z1 this 
regime corresponds to large tunnel couplings, 7 > Ao- 

We conclude by summarizing our qualitatively new 
finding: The proximity induced gap in SM-SC hybrid 
nanostructures is strongly suppressed in the intermedi- 
ate and strong tunnel coupling regimes whenever the SM 
layer thickness exceeds a characteristic crossover value 
determined by the band parameters of the SM. This leads 
to a strong coupling regime characterized by a small in- 
duced gap that is almost tunnel coupling independent, 
as actually observed experimentally. If the SC proxim- 
ity effect in the SM-SC hybrid structures engineered to 
host the elusive Majorana mode is indeed in the strong 
coupling regime, a critical reevaluation of the relevant 
theoretical predictions is urgently needed. 
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